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ABSTRACT 

We use two simulations performed within the Constrained Local UniversE Simulation 
(CLUES) project to study both the shape and radial alignment of (the dark matter compo- 
nent of) subhaloes; one of the simulations is a dark matter only model while the other run 
includes all the relevant gas physics and star formation recipes. We find that the involvement 
of gas physics does not have a statistically significant effect on either property - at least not 
for the most massive subhaloes considered in this study. However, we observe in both sim- 
ulations including and excluding gasdynamics a (pronounced) evolution of the dark matter 
shapes of subhaloes as well as of the radial alignment signal since infall time. Further, this 
evolution is different when positioned in the central and outer regions of the host halo today; 
while subhaloes tend to become more aspherical in the central 50% of their host's virial ra- 
dius, the radial alignment weakens in the central regime while strengthening in the outer parts. 
We confirm that this is due to tidal torquing and the fact that subhaloes at pericentre move too 
fast for the alignment signal to respond. 

Key words: methods: A^-body simulations - methods: numerical - galaxies: formation - 
galaxies: haloes 



1 INTRODUCTION 

Over the past decade there has been a tremendous increase in 
activity dedicated to satellite galaxies - the small, over dense, 
clumps of matter that inhabit dark matter haloes. Specifically, 
the study of satellite propertie s has been atta c ked from both 
a numerical perspect i ve (e .g. Gao et all 20041; Kravtsov et al 




2004 1: iGill et all 120041 120051: iDiemand et al.l 120071: [Kuhlen et al 
2003; IPereira et al.l 120081: IWarnick & Knebej 120061: ISales et al 
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Dolag et al.l 2009 ) as well as an 



observational one 



20091 
( e.g. [Xzzaro et alj 



2006 



Faltenbacher et alj 



20071: IKodosov et alj|2008t IWang et all 120081 : iBailin et alj|2008l : 
Chakrabarti & Blitj2009l:IWang et al.l2009h. The advent of "N ear- 
Field Cosmology" (cf. iFreeman & Bland-Hawthornl |2002|) has 
opened a new field dedicated to scrutinizing the (commonly ac- 
cepted) standard model of comology that is based on the exis - 
tence of dark matter and dark energy (e.g. iKomatsuet "alj |2009l) . 
Characteristics that were previously beyond modeling include, 
for example, the sh a pes, orientation, and alignment of subhaloes 
l Knebe et alj 120041: Izentner et all 120051; ILibeskind et alj 120051: 



Kuhlen et alj 120071 : iKnebe et alj 120081 : iFaltenbacher et al 
Pereira et al ] |2008h . Furthermore, while there is a clear 



2008; 



consensus 



that baryonic physics will have an effect on structure formation, 
especially on galactic and su b-galactic scales th a t are only now 
withi n reach of simulations dMaccio et alj 120061: IWeinberg et alj 
20081: lOkamoto et alj|2009l ; iTissera et alj|2009l : iDolag et 31.1120091 
Libeskind et al .1201 Oh . a rigorous quantification of these effects has 
not yet been fully completed. 

In this paper we use a set of two self-consistent cosmolog- 
ical simulations of the formation of the Local Group of galax- 
ie£| In one simulation we focus just on the gravitational in- 
teraction of dark matter whereas the second run incorporates 
baryonic physics as well. The simulations are studied with re- 
spects to their (dark matter) shapes as well as the so-called ra- 
dial alignment of subhaloes: this is the tendency of subhaloes to 
have their intrinsic shapes aligned towards the centre of their re- 
spective host halo. While such a phenomenon has been reported 
on a broad range of scales including massive clusters as well 
as galactic systems jHawlev & Peebles! 19751: Thompson! 19761: 
Pereira & K uhn 2005; Agustsson & Brainerdl2006l : IWang & White! 
20071 : Faltenbacher et alj |2008|) the literature still lacks the con- 
firmation of such a signal in gasdynamical simulations of cosmic 



1 See the "Constrained Local UniversE Simulations' (CLUES) project 
http : //www. clues-project . org for more details. 
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'"gas.i '"dm £ tyta 

DM: - 2.54 x10 s 0.15 5.29 X 10 7 

SPH: 4.42 X 10 4 2.1 X 10 5 0.15 5.29 X 10 7 



Table 1. Simulation parameters. From left to right, the initial mass per gas 
particle (in units of /i~'M ); mass of high resolution dark matter particles 
(in units of /i~'M G ); softening length (in units of /r'kpc); Number of high 
resolution DM particles. 



structure formation. Further, it has been argued that the shapes of 
(field) dark matter ha l oes can be modified by g a laxy formation (e.g . 



Katz & Gunnl Il99ll: iKazantzidis et alj 12004: iBailin et al. 




20051; 


Maccio et alj 120061: IWeinbere et alj 20081: Libeskind et al 
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eve n though the consequence may be quite small for dwarf galaxies 
(cf . | Diemand & Moorell2009l) and subhaloes (cf. Kazantzid is et alj 
1200% It therefore remains interesting to check whether or not there 
is a measurable effect for subhaloes and if the predicted (and ob- 
served) radial alignment signal persist when including such physics 
in the simulations. 

Please note that in this paper we investigate the dark matter 
shapes and the influence of baryonic physics on them; an in-depth 
study of the shapes of the baryonic component will be presented in 
a companion paper (Libeskind et al., in preparation). 



2 THE SIMULATIONS 

In this Section we describe the simulations used throughout this 
study and the methodology employed to identify host haloes and 
their substructure. 

2.1 Constrained Simulations of the Local Group 

In this paper, we use the s ame set of simulations as those published 
in iLibeskind etal] d2010h and refer the reader to that paper for a 
more exhaustive description and discussion of the various aspects 
of the simulations including the constraints that were imposed in 
order to reproduce a z = Local Group. As mentioned earlier, these 
simulations form part of the CLUES project. We briefly summarize 
the main properties of the simulation here, for clarity. 

We choose to run our simulations using standar d ACDM ini- 
tial c onditions, that assume a WMAP3 cosmology l lSpergel et alj 
|2007|) . i.e. n,„ = 0.24, Q b = 0.042, f2 A = 0.76. We use a normal- 
ization of cr 8 = 0.73 and an = 0.95 slope of the power spectrum . 
We used the PMTree-SPH MPI code GADGET2 l lSpringell [20051) 
to simulate the evolution of a cosmological box with side length 
of Lbox = 64/7~'Mpc. Within this box we identified (in a lower- 
resolution run) the position of a m odel local group that cl osely re- 
sembles the real Local Group (cf. ILibeskind et alj[2010h and re- 
simulated a 2/i~'Mpc region about the centre of this local group us- 
ing our full resolution equivalent to 4096 3 effective particles. Struc- 
tures that are still linear at z = are constrained to have various 
properties that are similar to objects in the local Universe, allowing 
us to simulate our host halos in an environment similar to the real 
one. 

We use the same set of initial conditions to run two simula- 
tions, one with dark matter only and another one with dark matter, 
gas dynamics, cooling, star formation and supernovae feedback. 
We shall call these two simualtions DM and SPH, respectively; 
the numerical parameters of our simulations are summarized in Ta- 
bled] 



For the gas dynamical SP H simulation, we follow the feed- 
back and star formation rules of ISpringel & Hemquistl d2003l) : the 
interstellar medium (ISM) is modeled as a two phase medium com- 
posed of hot ambient gas and cold gas clouds in pressure equilib- 
rium. The thermodynamic properties of the gas are computed in 
the presence of a uniform but evolving ultra-violet cosmic back- 
ground generated from QSOs and AGNs and switched on at z = 6 
Maardt & Madaulll996l) . Cooling rates are calculated from a mix- 
ture of a primordial plasma composition. No metal dependent cool- 
ing is assumed, although the gas is metal enriched due to super- 
novae explosions. Molecular cooling below 10 4 K is also ignored. 
Cold gas cloud formation by thermal instability, star formation, the 
evaporation of gas clouds, and the heating of ambient gas by super- 
nova driven winds are assumed to all occur simultaneously. 



2.2 The (Sub-)Halo Finding 

In order to identify halos and subhalos in our simulation we have 
ru n the MPI+OpenMP hybrid halo finder AHlQ described in detail 
in lKnollmann & Knebel J2009I) . AHF is an improvement of the MHF 
halo finder dGill et alj 2004 ), which locates local overdensities in 
an adaptively smoothed density field as prospective halo centres. 
The local potential minima are computed for each of these density 
peaks and the gravitationally bound particles are determined. Only 
peaks with at least 20 bound particles are considered as haloes and 
retained for further analysis. We would like to stress that our halo 
finding algorithm automatically identifies haloes, sub-haloes, sub- 
subhaloes, etc. For more details on the mode of operation and actual 
functionality we refer the r eader to the code description paper by 
iKnollmann & Knebel f200gh . 

For each halo, we compute the virial radius R^o, that is the 
radius r at which the density M{< r)/(4;rr 3 /3) drops below A v j,pb. 
Here pb is the cosmological background matter density. The thresh- 
old A v ir is computed using the spherical top-hat collapse model and 
is a function of both cosmological model and time. For the cosmol- 
ogy used in this paper, A v ; r = 355 at z = 0. Subhaloes are defined as 
haloes which (fully) lie within the virial region of a more massive 
halo, the so-called host halo. As subhaloes are embedded within 
the density of their respective host halo, their own density profile 
usually shows a characteristic upturn at a radius R, i, ZJhaio* where 
^haio would be their actual (virial) radius if they were found in iso- 
lation)^ We use this "truncation radius" R, as the outer edge of the 
subhalo (and refer to it as R sat ) and hence subhalo properties (i.e. 
mass, density profile, velocity dispersion, rotation curve) are calcu- 
lated using the gravitationally bound particles inside the truncation 
radius R,. For a host halo we calculate properties using the virial 
radius i?haio referred to as 2?h st- 

We construct merger trees by cross-correlating haloes in con- 
secutive simulation outputs. For this purpose, we use a tool that 
comes with the AHF package, called MergerTree that follows both 
host and subhalos identified at redshift z = backwards in time. At 
a previous snapshot, the code identifies as the direct progenitor the 



2 AMIGA halo finder, to be downloaded freely from 
http: //www.popia. ft .uam. es/AMIGA 

3 Please note that the actual density profile of subhaloes after the re- 
moval of the host's bac kground drops faster than for isolated haloes (e.g. 
IKazantzidis et alj2004l) ; only when measured against the main halo's back- 
ground will we find the characteristic upturn used here to define the trunca- 
tion radius r,. 
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shape study: radial alignment study: 

M sub > 2 X IO'/t'Mq M sub > 2 X 10 8 /i-'M o & b/a < 0.9 



DM: 


64 = 


17 + 47 


56 = 


16 + 40 


SPH: 


45 = 


20 + 25 


36 = 


18 + 18 



Table 2. The number of subhaloes in each model according to the various 
selection criteria. The splitting into two summands gives the number of sub- 
haloes closer (first number) and farther away (second number) than 50% of 
the host's virial radius. 



halo that shares the most particles with the present halo and is clos- 
es t in mass. Again, for m ore elaborate details we point the reader 
to lLibeskind et alj d2010h . 

These merger trees are in turn used to define the infall time 
of a subhalo, i.e. in our case the time it is fully embedded within 
the virial radius of its host: d sat +R S!lt < /fhost- The accuraacy of the 
determination of this infall time is sensitive to the time sampling of 
our outputs which are equally spaced in At a lOOMyrs for the DM 
model and At a 30Myrs for the SPH model, respectively. 



3 RESULTS 

Using our two simulations and their respective (sub-)halo analy- 
sis we study the shapes and the orientation of the shapes of sub- 
haloes with respect to their position within the host halo. We re- 
strict this analysis to the three most massive host haloes, namely 
the Milky Way, M31, and M33, and stack all subhaloes found in 
these hosts in all analysis and plots. In order to compare like with 
like, we make a mass cut, c omparing only the mos t massive of the 
z = subhaloes: following llibeskind et al. - who studied 

the same set of simulations with respect to the radial distribution 
of subhaloes - we consider only subhalos with masses larger than 
M sub > 2x 10 8 /r'M G , which roughly corresponds to subhaloes with 
more than 1000 particles. This cut was chosen in order to obtain the 
correct number of visible satellites orbiting about the Milky Way 
(~ 20). Further, for the analysis of the radial alignment in Sec. 13.21 
we wish to study satellites with well defined shapes and thus em- 
ploy the additional criterion b/a < 0.9 (with a > b > c being 
the eigenvalues of the m oment of iner tia te nsor) as su ggested (and 
succes sfully applied) in lPereira et al .1 (e.g. I2008T ) and lKnebe et al.1 
J2008I) . The respective numbers of subhaloes left for the analysis 
are summarized in Table [2] Note that we also split the sample ac- 
cording to whether the subhalo is within half the host's virial radius 
(4sat < 0.5 tfhost) or not (0.5 R hoit < d sat < R hoa ); the number of sub- 
halos falling into these two bins is also shown in Table [2] 



3.1 The Shapes of Subhaloes 

3.1.1 Determining the shape 

Before quantifying the shape and the radial alignment of shapes, we 
first must confirm that the subhaloes we identify, are in fact aspher- 
ical. Several options for calculating a subhalo's shape exist which 
we wish to investigate. To this extent we calculate the subhaloes' 
sphericity s = c/a, where a > b > c are the eigenvalues of the 
moment of inertia tensor of the subhalo. 

There are two commonly used definitions for this tensor. 
These are the (standard) moment of inertia tensor 
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Figure 1. The radial distribution of the medians of the sphericity .s = c/a for 
all particles (solid black) and the dark matter particles only (dotted black) 
in the SPH model. The thick lines correspond to a > b > c derived from the 
standard moment of inertia tensor while the thin lines utilized the reduced 
moment of inertia tensor. The error bars are the upper and lower quartile of 
the plotted medians. 



h.j / i ttln^jjiX jji (1) 
n 

and the reduced moment of inertia tensor 

$ = 2 m " x j -" I ' " ' ^ 

n 

where in both cases x,_„ and xp, are the ith and jth component of 
the nth particle coordinate and r„ its distance to the subhalo cen- 
tre; note that the reduced version puts more weight onto the central 
region; or ino ther words, the standard moment of intertia tensor 
weights each particle's contribution by r\ which is compensated 
for by the 1 /r\ factor in the reduced one. We calculate both Ijj and 
using all (bound) particles within the radius R sM of the sub- 
halo (or interior to some radius r < R sM when studying the radial 
shape profile). Therefore, our determination of the eigenvalues used 
for the shape determination is based upon a diagonolisation of the 
moment of intertia tensor based upon particles within a spherical 
sphere. This appr oach biases the shapes tow ards higher spherici- 
ties. As shown bv lBailin & Steinmeta d2005l) this can be corrected 
for, i.e. (c/a) = (c/a) S pJ erica] - We apply this correction to all shapes 
and refer to these corrected shapes as c/a throughout the paper. 
Note that using spherical shells rather than the ellipsoids does not 
affect the orientation of the principal axes and is hence irrelevant 
for the study of the radial alignment. 

In order to properly calculate subhalo shape we wish to de- 
termine how these two definitions for the moment of Inertia Ten- 
sors differ - not only between each other but also, for the SPH run, 
when considering all subhalo particles or just the dark matter com- 
ponent. In Fig.[TJwe show the median sphericity s = c/a, calculated 
using the two moment of inertia tensors, for the dark matter subha- 
los as well as for all particles in the SPH subhalos and the dark 
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matter component of SPH subhalofl The difference in the two def- 
initions is clearly seen: the reduced moment of inertia tensor (thin 
lines) is dominated by the inner region distribution well towards 
the edge of the subhalo; that is, the sphericity remains nearly con- 
stant at the peak value reached at approximately 0.1/? Kat . We further 
observe that for the SPH subhalos, the distributions for the total 
(solid back) and dark matter (dashed black) are practically indis- 
tinguishable, leading to the conclusion that baryons play a minor 
role in the shape determination of the subhalo - at least in the outer 
parts r/R m , > 0. 1. Despite the quality of our simulations we are 
still unable to resolve the very inner parts of these subhaloes. How- 
ever, we also observe (though not shown here) the drop towards 
str onger asphericity in the very central regions as previously found 
by iKuhlen et al.l J2007h in their dark matter only simulations. We 
additionally like to mention (also not explicitly shown here) that 
the drop in sphericity towards the subhalo edges (seen when utiliz- 
ing the standard moment of inertia tensor) is driven by subhaloes 
closer to the host's centre as these particular subhaloes feel stronger 
tidal forces. 

In Fig. Q] we showed the difference in the shape of subhaloes 
using either all particles or just the dark matter particles in the SPH 
run alone. We are therefore unable to draw any conclusions regard- 
ing the affect of baryons on an individual (subhalo's) shape, since to 
do this would require comparing the difference in shape between a 
given DM subhalo and the same subhalo modeled with SPH. How- 
ever Fig. [T] shows that it makes little (if any) difference if shapes 
are calculated using all or just the dark matter particles in the SPH 
model. 

However, one may rightfully raise the question whether or not 
our results are mere numerical artifacts or refer to a physical phe- 
nomenon. To shed some light on this issue we checked for two 
things: first, do we actually resolve the (inner parts of the) sub- 
haloes with sufficient baryonic particles; and second, are our results 
biased by the halo finder and the position of the subhalo within the 
host, respectively. 

In order to understand the lack of difference in the shapes 
when using either the total or dark matter only component in the 
SPH model, we need to examine the baryon fraction. We can con- 
firm that the median ratio of number of baryonic particles to the to- 
tal number of particles within a 0. l/f sa t is close to 35%. This ratio 
drops to about 20% at the subhalo's outer edge indicating that there 
is in fact a sufficient number of baryonic particles present through- 
out the subhaloes' volume. Note that this test confirms the presence 
of baryonic particles; their contribution to the shape/sphericitiy is 
further dependent on the mass (not number) fraction due to the 
weights nij in the definition both moment of inertia tensors, i.e. 
the actual baryonic mass fractions at the aforementioned radii are 
ft = 0.09 at 0. IR SM and f b — 0.055 at the outer subhalo edge, re- 
spectively. 

Additionally, our spherical overdensity halo finder AHF biases 
the shape determination owing to the fact that subhaloes closer to 
the host c entre will have sm aller radii (cf. also discussion in Sec- 
tion 3.2 in lKnebe et alj2008h . This is due to the embedding of sub- 
haloes within a sea of host particles such that objects closer to the 
central density peak of the host will be truncated earlier leading 
to a smaller radiues. To check whether or not our sample is af- 



fected by this effect we studied the relation between subhalo radius 
if S at and position d slli within the host. We do not find any signif- 
icant correlation for the subhalo sample under investigation with 
the Spearman rank coefficient of the two quantities being 0. 1 1 . To 
gauge the value of the Spearman rank coefficient we tested several 
pure random correlations consisting of the same number of parti- 
cles as subhaloes studied here and found the mean Spearman rank 
coefficient for these random samples to be 0. 10; hence our value of 
0.11 suggests no correlation]^ 

Lingering on the results obtained so far for a moment leads 
to the following conclusions: since the standard moment of iner- 
tia tensor is unbiased by the central region and thus more sensitive 
to (radial) shape variations, we will use this measure rather than 
the reduced moment of inertia tensor throughout this work. Or put 
differently, we are interested in the shape variation as a function 
of radius r/R silt where the shape at a particular radius should be 
calculted utilizing a shell of thickness [r/R sdt - e, r/R sat + e], How- 
ever, despite the quality of our simulations there are still not enough 
particles within such shells for the subhaloes and hence we prefer 
to employ the standard moment of inertia tensor that actually puts 
more weight onto the particles in the region close to r/R sa . Further, 
as baryons have no effect on the subhalo shape determination, we 
use all subhalo particles for the shape determination in both our 
models (i.e. DM and SPH). 

We wish to remind the reader that up to now none of the results 
suggests that the baryons do not have an influence on the shape as 
this requires a direct comparison of the DM vs. SPH model (see 
Subsection l3.1.3l i. 



3.1.2 Radial profile of host sphericity 

Even though we primarily address subhalo shape, we neverthe- 
less like to also present the (mean) radial sphericity profile for 
the three host haloes under investigation, too; this allows us to 
get a better feeling and gauge for the differences between sub- 
haloes and host haloes, respectively. This subject has be en the tar 



get of many previous investigations by other groups ( e.g. Dubinski 
19941 : iGustafsson et all 120061 ; iDebattista etal] 120081 ; lAbadi et al 
200S jTissera et alj|20090 with the conclusion that the inclusion of 



baryons drive the dark matter halo to become more spherical. 

Our results (based upon the standard moment of inertia tensor) 
can be viewed in Fig. [2] where we show the sphericity for the DM 
model (dot-dashed line) alongside the two profiles, i.e. DM only 
(dashed line) and total matter (solid line), for the SPH model. We 
clearly see differences in the sphericities in the SPH model below 
0.1i?h ost in the sense that the dark matter component alone appears 
more spherical than the total ma tter. This transition at 0.1i?h ost coin- 
cides with the value reported bv lBailin et all j2005>) who found that 
below this scale the inner halo appears to be correlated with the 
disk in their simulations. However, we also note that the DM only 
simulation shows a rather strikingly different radial sphericity pro- 
file indicative of a substantial influence of the baryonic material on 
the (central) shape: even at 10% of the host's virial radius there is 
a difference of approximately 0.12 with the DM model being more 



4 Note that we are not utilizing the shape of the baryonic component in 
this particular study; we are mainly interested in the (possible) differences 
in the dark matter component arising from the inclusion of gas physics into 
the cosmological simulation. 



The Spearman rank coefficient is a non-parametric measure of correla- 
tion: it assesses how well an arbitrary monotonic function describes the 
relationship between two variables, without making any other assump- 
tions about the particular nature of the relationship between the variables 
iKendall & Gibbonslll990l) . 
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Figure 2. The mean radial distribution of the sphericity of our three host 
haloes in the DM (dot-dashed line) and SPH model (solid and dashed lines). 
The measurements are based upon the standard moment of inertia tensor 
and for the SPH model upon the total (solid line) and DM only component 
(dashed line). 

aspherical. The situation actually reverses when moving even fur- 
ther in and below 0.05/?host the dark matter component in the DM 
model now appears more spherical than in the SPH simulation; a 
clear manifestation of the baryonic mass concentration in the cen- 
tre and the disk structures seen in all three host haloes, respectively. 
Another important feature to note is that for the host haloes we do 
see a difference between the shapes of the DM and total matter 
component as opposed to the subhaloes (cf. Fig. [I}: they developed 
stronger and more concentrated baryon accumulations in the centre 
to actually affect the dark matter. This is confirmed by calculating 
again (though not presented) the radial dependence of the ratio of 
number of baryonic to total particles for the hosts: this ratio is ac- 
tually approximately 40% at /?host> about 75% at 0. lR\ msi and even 
higher in the very central regions. Again, we prefer to measure the 
influence in "number of particles" as the usage of more particles 
entails a better sampling of (the radial dependence of) shapes. For 
completeness, the baryonic mass fractions at aforementioned host 
radii are f b = 0.2 at 0. llfhost and f b = 0. 1 at /fhost. respectively, with 
that fraction increasing to more than 0.8 in the central regions. 

Throughout the rest of this paper we undertake the same com- 
parison for the subhaloes, i.e. DM vs. SPH, though restricting the 
analysis of the SPH model to the total matter component as we do 
not find differences in the total and dark matter shapes (cf. Fig.[T}. 

3.1.3 Radial profile of subhalo sphericity 

Having established how to calculate the shapes (standard moment 
of inertia tensor) and which component in the SPH run to use (total 
matter content) we now turn to a comparison of the shapes of the 
DM and SPH subhaloes. In Fig.[3]we present the radial distribution 
of the median (alongside ± 25%) values of the sphericity for the 
DM (black) and SPH (red) model both at redshift z = (solid lines) 
and at the respective infall time of the satellites (dashed lines). The 
most striking feature of this plot is the lack of any (statistically 
significant) difference between the DM and SPH runs, irrespective 
of the redshift. This implies that while baryon physics may affect 
the shapes of host haloes (e.g. lBailin et al . 2005 1 and Section [3X2l , 
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Figure 3. The radial distribution of the medians of the sphericity s = c/a 
based upon a > b > c derived from the standard moment of inertia tensor 
for the DM and SPH model, respectively. The black lines correspond to the 
dark matter only simulation while the red lines are measures including all 
particles (dm, gas, and stars) in the SPH run. The solid lines refer to redshift 
z = while the dashed lines are based upon the subhaloes' sphericities as 
measured at their respective infall times. 

there appears to be little influence on subhaloes. Recall that we 
showed in Fig. [2] that we do in fact observe an influence of the 
ba ryons on the shape of t he host haloes (in agreeement with Fig. 2 
in iLibeskind etal] fcOld) where we showed the influence of the 
baryonic component on the density profile inwards from approx. 

O.lflhost). 

However, we find that the shape evolves from the time a sub- 
halo fell into the host and the present time: subhalos in the central 
parts become rounder while those in the outer parts become more 
apsherical . The latter is readily explained by tidal forces to be in- 
vestigated in more detail in the following subsection. 

3.1.4 Dependence of sphericity on subhalo position within host 

To study the importance of tidal forces on shape, we split our sam- 
ple into subhaloes closer and farther away than 50% of the host's 
virial radius and present the sphericity profiles for each subset in 
Fig. [4 ] Note that th is division is motivated by the observations of 
IPereira et aljj2008l) which indicate that the peak of the radial align- 
ment as a function of distance from the host's centre lies at approx- 
imately half the host's virial radius (cf. their Fig.4). The upper and 
lower plot refer to subhaloes closer and farther away than 0.5R<n OSt 
at redshift z = 0. Note that this division is based upon the dis- 
tance d sM at today's time; at infall time every satellite obviously has 
d»a ~ ^host- We find that the sphericity change is primarily driven 
by objects closer to their respective hosts as expected if tidal forces 
were responsible for this result. We also observe a tendency for the 
closest SPH subhaloes to become more spherical than their DM 
counterparts in the central regions while there is a marginal oppo- 
site trend fo r the outer subhaloes in a greement with the findings of, 
for instance jKazantzidis et al.l J2004I) . 

To better quantify the sphericity change (as measured for all 
our subhaoes at their respective radii R sat ) we plot in Fig. [5] the 
ratio s z=0 /s z=Zi against the position d siit of the subhalo within the 
host, where s :=0 is the sphericity as measured today and s z=zi at 
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Figure 4. The same as Fig.[3]but this time dividing the subhaloes into two 
populations: those that end up at z = closer to the centre than 50% of their 
respective host halo's radius (upper panel) and the ones between 50-100% 
of the host's radius (lower panel). Note that this division is based upon 
the present day position d sM of the subhalo; at infall time every satellite 
obviously has d sal a Rhost- 

infall time. This ratio is primarily below unity indicating a drop in 
sphericity. Further, subhaloes ending up closer to their host's centre 
show a (marginally) stronger drop in agreement with the picture of 
tidal origin. This correlation, (i.e. stronger sphericity evolution for 
subhaloes closer to their host) is confirmed by the Spearman rank 
coefficient taking a value of 0.46 for both the combined DM and 
SPH model and a linear fit to the combined data giving a slope of 
0.35 (indicated by the solid line). 



3.2 The Radial Alignment 

As we have just seen, subhaloes/satellites themselves are - just 
like field ha loes - asp herical systems as previously shown by 
I Knebe et al .1 (e.g. 1200 8h . Furthermore, it has also been reported 
in the literature that the orientation of their shapes is not random 
and the ir major axis preferentially points towards their host cen- 
tre (e.g. lKuhlen et alj|2007l . lFaltenbacher et alj|2008l ; IPereira et al.l 



* DM 




best fit linear slope 0.35 + 

i i i i i i i i i i I i i i L_ 

0.0 0.2 0.4 0.6 0.8 1.0 

^sat/^-host 



Figure 5. The ratio of sphericities .v : at redshift z. = and s Zi at the 
respective infall redshift Zi of the subhalo as a function of position J sat of 
the subhalo at redshift z = within the host's virial radius R n ost ■ The solid 
line indicates the best fit linear relation for the combined sample (i.e. DM 
and SPH together). 

120081 ; iKnebe et alj|2008lfl) . This "radial alignment" is quantified 
by the (cosine of the) angle between the position vector of the sub- 
halo cUat in the rest frame of its host halo and the major axis a sat of 
the actual subhalo (again, as determined via the eigenvectors of the 
standard moment of inertia tensor). This quantitiy 

cos (/> = d sal ■ a sal (3) 

will be the target of our study in the following subsection, in order 
to be able to properly measure th e radial alignment, and in com- 
plianc e with previous studies (cf. IPereira et alj |2008; Kneb e et"al] 
l2008h . we wish to further restrict ourselves to those subhaloes with 
well defined major axes, pruning our sample such that all subhaloes 
have an intermediate to long axis ratio of b/a < 0.9. Please refer to 
Table|2]for the number of subhaloes in this restricted sample. 



3.2.1 Radial profile of radial alignment 

The radial alignment (and its evolution since infall time) is shown 
in Fig. [6] where the radial distribution of the medians of |cos^| as 
defined by Eq.((3j is plotted. We note that while the signal appears 
to be strong (as expected) there is no noticeable difference between 
the DM and SPH model. This comes as no surprise as we previ- 
ously mentioned that there is virtually no difference between the 
shapes of subhaloes in the DM and the SPH model (cf. Section [3~TT l. 
That said, we note the trend for the signal to strengthen over time - 
akin to the observation of the evolution of sphericities. This is ex- 
pected if the signal ste ms from evolutionary r ather than env i romen - 
tal effects as show n by iKuhlen et al] d2007f) ; IPereira et ail ( 120081) : 
IKnebe et al.l d2008l) . We further find that the signal weakens when 
moving the point of reference for the sub halo shape inside, just as 
already indicated bv lKuhlen et alj J2007l) . However, this cannot be 
related to the fact that the subhaloes become more spherical in the 
inside (cf. Fig. [3) as (though not explicitly shown here) the radial 
aligment follows a similar radial trend throughout the subhalo when 
using the reduced moment of inertia tensor. 
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Figure 6. The radial distribution of the medians of the radial alignment 
based upon a > b > c derived from the standard moment of inertia tensor. 
The black lines correspond to the dark matter only simulation while the red 
lines are measures including all particles (dm, gas, and stars) in the SPH 
run. The solid lines refer to redshift z = while the dashed lines are based 
upon the subhaloe's sphericities at their respective infall times. 



3.2.2 Dependence of radial alignment on subhalo position 
within host 

As before, we divide our subhaloes into two distinct samples ac- 
cording to position in the halo: objects closer to the centre than 50% 
of their host's virial radius and satellites between 50 and 100% of 
the host's radius. This division is shown in Fig. UJ where the dis- 
tribution of the radial alignment for the inner (outer) sample can 
be viewed in the upper (lower) panel. We find that the signal pri- 
marily evolves for those subhaloes that remain in the outskirts of 
the halo: subhaloes that end up closest to the host's centre at red- 
shift z = experience a mild weakening of the radial alignment - 
in fact the evolution is also consistent with no evolution due to the 
large scatter. At first this may appear surprising as we previously 
showed that the shapes become more aspherical over time primarily 
for those objects ending up closer to the host. One may thus expect 
a strengthening of the radial alignment. But as we show later, the 
scenario is a bit more complicated and is related to the fact that sub- 
haloes m ove faster at pericent re. However, our finding s are consis- 
tent with lKuhlen et all J2007I) and lPereira et al.l d2008h who found 
a considerably stronger signal for satellites between [0.5, 1.0]^ n ost 
than for objects closer than this. These authors explain this by as- 
serting that the subhalo is moving too fast at pericent re for tidal 
torquing to be effective (cf. Fig. 8 in lPereira et alj|2008l) . Below, in 
subsection l3.3l we provide direct evidence for this scenario. 

We conclude this subsection with an investigation of the cor- 
relation between the radial alignment and subhalo position on an 
individual object basis. To this extent, we present in Fig. [8]the ra- 
dial alignment cos <p = d sal ■ a sal for each subhalo vs. its position 
c4at as measured today (upper panel) and the ratio between d sal ■ a sa t 
as measured today and at infall time (lower panel). We observe 
that the radial alignment at z = shows a (marginal) tendency to 
be stronger for objects in the outskirts of the host halo - as con- 
firmed by reading off the corresponding numbers for the medians 
at r/R sM = 1 in the upper and lower panel of Fig. [7] The lower panel 
of Fig.[8]further confirms a strengthening of the signal over time - 



1.0 




r/R, 



Figure 7. The same as Fig.|6]but this time dividing the subhaloes into two 
populations: those that end up at z = closer to the centre than 50% of their 
respective host halo's radius (upper panel) and the ones inbetween 50-100% 
of the host's radius (lower panel). 

at least for the subhaloes residing in the outer parts of their host 
halo at redshift z = as noted earlier. This correlation between 
evolution of radial alignment and subhalo position in the host is 
again confirmed by the rather large value of 0.74 for the Spearman 
rank coefficient as well as the best fit linear slope of 0.57. Please 
note that the lower panel also indicates a (mild) weakening of the 
radial alignment for objects ending up close to the host in agree- 
ment with the picture that the signal did not have enough time to 
adjust itself for such "speeding" subhaloes. We shall have a closer 
look at exactly this phenomenon in the subsequent Subsection [33] 

3.3 Correlations with Subhalo Velocity 

One explanation for the stronger radial alignment in the outer parts 
is tidal torquing and the fact that subhaloes are simply moving too 
fast at pericentre f or the radial alignment to become effective (e.g. 
IPereira et alj|2008h . To confirm this picture we examine the corre- 
lation of the long axis a sat with the subhalo's velocity v sat in the host 
halo's rest frame: 
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Figure 8. The radial alignment as given by Eq.{3} as a function of the sub- 
halo's position rf sat within its host halo. The upper panel shows the relation 
for redshift z = while the lower panel shows the ratio between the value 
at redshift z = and the respective infall redshift of the subhalo. 



Figure 9. The correlation between the major axis a sal of the subhalo (as 
measured at its radius K sat ) and the velocity v sat in the rest frame of the 
respective host halo as a function of position d wjX of the subhalo within its 
host. The upper panel shows the relation for redshift z = while the lower 
panel shows the ratio between the value at redshift z = and the respective 
infall redshift of the subhalo. 



cos (f> = v sat • a sat . (4) 

Fig.[9]shows the modulus of cos <p as a function of the distance 
d slit to the host centre both at redshift z = (upper panel) and the 
ratio between v sat • a sat as measured today and at infall time (lower 
panel). We note that subhaloes ending up in the central regions of 
their hosts have their shapes (as measured by the major axis a sat 
at R sM ) aligned with the velocity vector v sat . This explains why we 
saw a weakening of the radial alignment signal in Section [X2l the 
velocity vector is tangential to the orbit and hence not pointing to- 
wards the host. However, the shape appears to be elongated along 
the flight path (which we have confirmed by looking at d sat ■ v sal - 
though not explicitly presented here - where we found a (mild) cor- 
relation at infall time) and hence the correlation of both quantities. 
This obse rvation is naturally i n agreement with the previous expla- 
nation of IPereira et alj ( 120081) : the subhalo does not have enough 
time to re-arrange its shape at pericentre passage; it nevertheless 
gets stretched/tidally torqued when passing close to the host cen- 



tre. However, it appears that this signal is already present at infall 
time, i.e. even the infall shape and infall velocity are more strongly 
correlated for those subhaloes ending up closer to the host which is 
confirmed by the rather weak evolution of v sat • a sat as seen in the 
lower panel of Fig.[9] There we further find only a mild dependence 
of this evolution with position d sal within the host as confirmed by 
both the best fit linear slope of -0.17 and the Spearman rank coef- 
ficient of -0.28 for the combined sample. 

We like to close with a cautionary remark: from both Fig. [8] 
and Fig. [9] one may be inclined to infer that the (infall) velocity of 
subhaloes should be randomly distributed as there appears to be no 
correlation between shape (as measured by a sa( ) and v sal whereas 
the shape is correlated with the position vector d sat for the sub- 
haloes ending up in the outer parts of their respective host halo. 
However, one needs to bear in mind that subhaloes on radial or- 
bits (i.e. v sat aligned with d sat ) preferentially end up in the central 
regions; subhaloes lingering in the outer parts are favourably on 
more circular orbits. This picture is confirmed by actually plotting 
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v sat • d sal which we omitted here: those subhaloes who had both 
these vectors aligned at infall time (i.e. radial orbits) actually end 
up closer to the centre with smaller rf sa t/^hosi; the ones remaining 
in the outskirts show a weaker alignment at infall time. 



4 SUMMARY AND CONCLUSIONS 

In this paper, we use two simulations performed within the CLUES 
project to study both the shape and radial alignment of subhaloes; 
one of the simulations is a dark matter only model while the other 
run includes all the relevant gas physics and star formation recipes. 

We first set out to determine the most accurate method to mea- 
sure the shape of subhaloes. To this end we compared two similar 
methods: the standard and reduced moment of inertia tensor. We 
found that due to the 1 jr 2 -weighing, the reduced moment of inertia 
tensor gives a stronger bias towards the central shape of subhaloes. 
Since in this work we are interested in the effects of tidal forces that 
primarily act in the outer parts, we decided to employ the standard 
moment of inertia tensor. We further investigated the dependence of 
measuring the sphericity of SPH subhaloes on the individual com- 
ponents, e.g. solely the dark matter particles or all particles includ- 
ing gas and stars. Despite the presence of baryonic particles in our 
SPH subhaloes, we are unable to find any significant differences in 
the sphericities of subhaloes when calculated using either all the 
particles or just the dark matter ones. We thus decided to use all 
particles when determining the shape of SPH subhalos. This does 
not automatically imply that baryons have no effect on shape just 
that in an SPH subhalo the baryonic and dark matter components 
follow the same radial shape profile. 

Our first result is that - when comparing the sphericities of 
subhaloes in the DM and the SPH model - the distribution of 
shapes is not affected by including baryons. We further note that 
we also investigated the triaxiality in the same way as the spheric- 
ity and were led to the same conclusions; therefore, we use the 
terms sphericitity and shape interchangeably. However, there is (for 
both models) an evolution in time: subhaloes become more aspher- 
ical over time with this evolution being driven primarily by objects 
closer to the centre of the host, as expected if tidal forces were re- 
sponsible for alterations to the shape. 

We further find that the radial alignment is not affected by gas 
physics, which should come as no surprise since we did not find an 
influence of gas physics on the shapes of subhalos. We thus confirm 
the picture of an evolutionary (rather than environmental) origin of 
the signal: there is considerable evolution of the radial alignment 
since infall of a subhalo with the signal being weaker for objects 
close to the centre. This weakening can be explained by the fact that 
the subhaloes move too fast at pericentre for the alignment signal 
to adjust. This is confirmed by a correlation between the shape of 
the subhalo and its velocity vector, the tangent to the actual orbit 
not pointing towards the host centre. 

The most remarkable outcome of this study though remains 
that the inclusion of gas physics has no affect on the (dark mat- 
ter) shapes (and hence radial alignment) of subhaloes in cosmolog- 
ical simulations of structure formation. While we do find an influ- 
ence in our suite of host haloes along the lines reported by other 

I : II 1 I II ^1 \ 1 

authors (e.g. Dubinski 1994; Bailin et al. 2005: Gustafsson et al. 

l2006l : lDebattista et alj2008l : lAbadi et al.l2009l : lTTs'sera et al.l200"9l) , 

there appears to be no impact on subhaloes. Why are subhaloes 
different to host haloes in that respect? After all, the inclusion of 
baryons should affect the dark matter, too: lBlumenthal etak 
showed that dissipative baryons will lead directly to the adiabatic 



contraction of haloes (see also lGnedin et alj[2004l) and should thus 
be a critical ingredient in determining halo and subhalo properties. 
However, here we find that while it most certainly affects the host 
haloes (cf. Fig. |2j it appears unimportant for the subhaloes. We 
though need to note that this paper looks primarily at low-mass 
subhaloes (Milky Way dwarf galaxies), and hence one could imag- 
ine the result changing significantly on cluster scales, where the 
substructure now corresponds to massive galaxies, where the dis- 
tribution and affect of baryons might be quite different. 

Although we are confident that our results are not driven by 
numerical artifacts (see Section [3.1.1t . we note that the fraction of 
baryonic particles in the host haloes is substantially larger than in 
our subhaloes: 20-30% of the particles in a given subhalo are bary- 
onic particles. This ratio increases with decreasing subhalo radius: 
the baryon number fraction is approximately 40% at i^ost, ~75% 
at 0. l.Khost and even higher in the very central regions. Note that 
we chose to adhere to "number fractions" as a a greater number of 
particles will allow for a better sampling of the (radial dependence 
of) (sub-)halo shape - irrespective of their mass; the correspond- 
ing baryonic mass fractions can be found in Sec tion [3. 1.1 l and Sec- 
tion [37L2] respectively. Given this difference in host and subhalo 
fractions we like to caution that even higher resolution may be re- 
quired to confirm the results presented here. We further checked 
whether the baryons in the subhaloes are primarily in the form 
of gas or stars. In that regards we found that the gas mass frac- 
tion is substantially lower than the stellar mass fraction. In fact, 
the median gas to stellar mass fraction is of order 5%, i.e. 95% of 
the baryons are stars. Under the assumption that those stars have 
formed prior to the infall of the subhalo into its host and the fact 
that star particles are effectively behaving like dark matter particles 
we may draw the conclusion that our results are not as suprising as 
one may think: even though the stars have a different formation his- 
tory, i.e. they are born out of gas which had the option to cool and 
sink to the centre of the subhaloes' potential wells, the influence 
of the baryonic/stellar content of subhaloes will be limited as ob- 
served in this paper. However, the study presented here focused on 
the shape of the dark matter only component. To further investigate 
the influence of the baryonic component we will in an upcoming 
paper, study the shapes of the baryonic component (Libeskind et 
al., in preparation). 
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